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Abstract 



In this paper we consider a discrete scale invariant Markov process {X{t),t G R"*"} 
with scale ^ > 1. We consider to have some fix number of observations in every scale, 
say T, and to get our samples at discrete points a'^, k ^ W, where a is obtained by the 
equality I = and W = {0, 1, . . .}. So we provide a discrete time scale invariant Markov 
(DT-SIM) process X{-) with parameter space {a^,k £ W}. We present some properties 
of such a DT-SIM process and we show that the covariance function is characterized by 
the values of (1), i?f (0), j = 0, 1, . . . , T - 1}, where Rf{k) is the covariance func- 
tion of jth and (j -I- fc)th observations of the process. We also define the corresponding 
T-dimcnsional self-similar Markov process and characterize its covariance matrix. 
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1 Introduction 

The notion of scale invariance or self-similarity is used as a fundamental property to handle 
and interpret many natural phenomena, like textures in geophysics, turbulence of fluids, data 
of network traffic and image processing, etc [1]. The idea is that a function is scale invariant 
if it is identical to any of its rescaled functions, up to some suitable renormalization of its 
amplitude. 

Discrete scale invariance (DSI), a property which requires invariance by dilation for certain 
preferred scaling factors, was pioneered by Sornette [10]. A random process has DSI if its 
finite dimensional probability distributions are globally invariant under the action of scaling 
operators of a fixed ratio [13]. It is known that DSI leads to log-periodic corrections to scaling. 
Log-periodic oscillations have been used to predict price trends, turbulent time series and 
crashes on financial markets [11]. Flandrin et. al. have studied the property of DSI and its 
relation to cyclostationarity by means of the Lamperti transformation [2], [5]. Discrete time 
linear systems that possess scale invariance properties even in the presence of continuous 
dilation were proposed by Zhao and Rao [T2]. DSI time series in real life usually have some 
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scale which is not necessarily integer. Thus we consider some special sampling scheme to 
obtain a discrete time DSI process, sampled from a continuous time DSI one. 

Markov processes have been the center of extensive research activities and wide sense 
Markov processes are studied before the general theory. In some texts, these processes are 
defined in the case of transition probabilities of a Markov process. Various classes of wide 
sense Markov processes are, like jump processes, diffusion processes and processes with a 
discrete interference of chance [6]. Measurements on a network have shown that LAN traf- 
fic behaves as a self-similar process. Through experiments, it has been shown that certain 
Markov processes can behaves as self-similar processes over a limited but large time do- 
main [9j. By applying both the Markov property and the self-similarity, we have self-similar 
Markov processes. These processes are involved in various parts of probability theory, such 
as branching processes and fragmentation theory. 

We consider a DSI process with some scale / > 1, and we get our samples at points a^, 
where k € W, / = , W = {0, 1, . . .} and T is the number of samples in each scale. By 
such sampling we provide a discrete time scale invariant Markov (DT-SIM) process in the 
wide sense which is a discrete time scale invariant process and is Markov in the wide sense 
too. We find the properties of such a DT-SIM process and also the covariance function is 
characterized. 

The paper is organized as follows. In section 2, we present definitions and some pre- 
liminary properties of Markov processes in the wide sense, DSI and self-similar processes 
by special operators. We define discrete time self-similar and discrete time scale invariant 
processes too. In section 3 we obtain a closed formula for the covariance function of the 
DT-SIM process. We also define and characterize the covariance matrix of corresponding 
T-dimensional self-similar Markov process in section 4. 

2 Theoretical framework 

In this section, by using renormalized dilation operator, we introduce self-similar and wide 
sense self-similar processes and define the concepts of discrete time self-similar and discrete 
time scale invariant processes. We also present some characterizations of Markov processes 
in the wide sense. 

2.1 Stationary and self-similar processes 

Definition 2.1 Given r E R, the shift operator Sr operates on process {Y{t),t G R} ac- 
cording to 




(2.1) 



A process {Y{t),t G R} is said to be stationary, if for any t, r € R 



{{SrY){t)}^{Y{t)} 



(2.2) 



where = is the equality of all finite- dimensional distributions. 
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If (2.2) holds for some r G R, the process is said to be periodically correlated. The smallest 
of such T is called period of the process. 

Definition 2.2 Given some numbers H > and A > 0, the renormalized dilation operator 
T^H,\ operates on process G R^} according to 

{VH,xX){t) ■= A-^X(At). (2.3) 

A process {X{t),t £ R^} is said to be self-similar of index H , if for any A > 

{{VH,xX)it)} ^ {X{t)}. (2.4) 

The process is said to be DSI of index H and scaling factor Aq > or (H,Ao)-DSI, if 
(2.4) holds for X = Xq. 

Definition 2.3 A process {X(k),k G T} is called discrete time self-similar process with 
parameter space T , any subset of distinct points of real line, if for any ki,k2 (zT 

{X{k2)} ^ (|)^{X(fci)}. (2.5) 

The process X{-) is called discrete time scale invariance with scale L > and parameter 
space T , if for any ki,k2 = Lki G T, (2.5) holds. 

Remark 2.1 // the process {X{t), t G R^} is DSI with scale I > 1, then for any fixed s > 0, 
X{-) with parameter space T = {I'^s^k G Z} zs a discrete time self-similar process. If the 
process {X{t),t G R^} is DSI with scale I = , for some T G N and a > 1, then for any 
fixed s > Q, X[-) with parameter space T = {a^s; k G Z} is a discrete time scale invariant 
process, with the same scale I = . 

Remark 2.2 // the process {X{t),t G R^} is DSI with scale I = for fixed T G N and a > 
1, then by sampling of the process at points a^, k G W where W = {0, 1, . . .}, we have X{-) 
as a discrete time scale invariant process with parameter space T = {a^; k G W} with scale 
I = . If we consider sampling of X{-) at points a"^"'"'^,n G "W, for fixed k = 0,1, . . . ,T — 1, 
then X{-) is a discrete time self-similar process with parameter space t = {a"'^+'=;n G W}. 

Definition 2.4 A random process {X{t),t G R^} is said to be wide sense self-similar with 
index H, for some H > if the following properties are satisfied for each a > 

(i) E[X'^{t)] < oc, 

(ii) E[X{at)] = a^E[X{t)], 

{Hi) E[X{ati)X{at2)] = a'^" E[X{ti)X{t2)]. 

This process is called wide sense discrete scale invariance of index H and scaling factor 
uq > 0, if the above conditions hold for some a = oq. 
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Definition 2.5 A random process {X{t),t G R^} is said to be discrete time self-similar 
(DT-SS) in the wide sense with parameter space T = {A^s, n £ W, A > 1} for any fixed s > 
and index H > 0, if the followings are satisfied for every n, ni, n2 G W 



(i) ^[X2(A"s)] < oo, 

{ii) £:[X(A"+"is)] = A"-f^^[X(A"is)], 

(m) ^[X(A"+"is)X(A"+"25)] = A2"-f^£;[X(A"is)X(A"2s)]. 



// the above conditions hold for fixed n, then the process is called discrete time scale 
invariance (DT-SI) in the wide sense with scale / = A". 

Through this paper we are dealt with wide sense self-similar and wide sense scale invariant 
process, and for simplicity we omit the term "in the wide sense" hereafter. 

2.2 Discrete time Markov processes 

Let {X{n),n G Z} be a second order process of centered random variables, E[X{n)] = and 
i?[|X(n)p] < cx), n G Z. Following Doob [1], a real valued second order process {X{n),n G Z} 
is Markov in the wide sense if, whenever ti < ... < tn, 



is satisfied with probability 1, where E stands for the linear projection (minimum variance 
estimator). If the process is Gaussian, then E is a version of the conditional expectation. 
The following facts on covariance of Markov processes in the wide sense, are essentially due 
to Doob !4l. Let 



where R{ni,n2) ■= E[X(ni)X{n2)] is the covariance function of X(-). Then {X{n),n G Z} 
is Markov if and only if g satisfies the functional equation 



E[X{tn)\X{ti), . . . = E[X{tn)\X{tn-l)] 




R{ni,ni) > 
R{ni, ni) = 



Q{ni,n2) = Q{ni,n)g{n,n2) 



ni ^ n ^ 722 



which is the same as 



R{ni,n)R{n, 112) 



R{n, n)i?(ni, 71-2) 



ni ^ n ^ n2. 



It follows that 

R{ni,n2) 
for some functions G and K [7j. 



Gini)K{n2) 



ni ^ 71-2 
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Borisov [3j completed the circle even for continuous time processes, namely, let R{ti,t2) 
be some function defined on T x T and suppose that R{ti,t2) 7^ everywhere on T x T, 
where T is an interval. Then for R{ti,t2) to be the covariance function of a Gaussian Markov 
process with time space T it is necessary and sufficient that 

Riti,t2) = G{mm{ti,t2))K{max{ti,t2)) (2.6) 

where G and K are defined uniquely up to a constant multiple and the ratio G/K is a positive 
nondecreasing function on T. 

It should be noted that the Borisov result on Gaussian Markov processes can be easily 
derived in the discrete case for second order Markov processes in the wide sense, by using 
theorem 8.1 of Doob [1]. 



3 Covariance structure 

In this section we study the structure of the covariance function of discrete time scale invariant 
Markov process and denote it by DT-SIM. Let {X{t),t E R-^} be a mean zero DT-SIM 
process with scale I. If / < 1, we reduce the time scale, so that / be greater than 1 in 
the new scale. Now we consider to have T samples in each scale. Our sampling scheme 
is to get samples at points 1, a, a^, . . . , a^^^, , a^~^^, . . ., where a is obtained by equality 
I = . Thus the process under study will be {X(a"),n G W} which is DT-SIM with scale 
I = c? . We show that the covariance function R^i ) of DT-SIM is specified by the values of 
{Rj^{l), Rj^ {0),j = 0, 1, . . . , T — 1}. We present a closed formula for the covariance function 
of DT-SIM process in theorem 3.1. The covariance function of DT-SIM process under some 
conditions is characterized in theorem 3.3. 

Theorem 3.1 Let {X(a"),n G W} be a DT-SIM process with scale I = , a > 1, T e N , 

then covariance function 

R^{t) = i?x(a"+^ a") := S[X(a"+^)X(a")] (3.1) 
where r G W, n = 0, 1,...,T— 1 and R^{t) is of the form, 

R^{kT + v) = [/i(a^-i)]^/i(a^+"-i)[/i(a"-i)]-^i?^(0) (3.2) 
where k G W, v = 0,1, . . . ,T — 1 

r r 

h{a') = Yl h{a^) = JJ Rf{l)/Rf{0), r G W (3.3) 
i=o i=o 

and h{a~^) = l.D 
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As {X{a^),n G W} is DT-SI with scale , this theorem fully characterize the covariance 
function of the DT-SIM process. Before proceeding to the proof of the theorem we present 
some equalities about the covariance function of the process which is needed in the proof of 
the theorem. 



Prom the Markov property (2.6), for a > 1, -R^(r) satisfies 

(r) = G(a")ii:(a"+^) r G W. 
By substituting r = in the above relation we have that 



(0) = Gia^Kia^) =^ G{a^) = 



Therefore 



" K{a'^) 



r G W 



i^(a"). 



Now it follows that 
Thus 



«1) 



n-l 



K{a'^) = K{1) Yl h{a^) 

j=0 



where h{a^) = Rf{l)/Rf{0). Hence forn = 0, 1, . . . , T - 1, keW 

fcT+n-l 

^(^feT+n) ^ ^^^^ -Q ^(^jy 
j=0 

As X{-) is DT-SI with parameter space {a*', k G W} and scale by (3.1) 

R?+.{i) _ i?f (1) 



R^^,{0) i?f(0)' 



i G W. 



Therefore 



fcT+n-l 

n M 

j=0 



a-' 



n -=0 ^(«^' 



n k 



n > 
n = 



(3.4) 



(3.5) 



(3.6) 



(3.7) 
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Thus using (3.3) and the convention h{a ^) = 1 we have 

fcT+n-l 



Yl h{a^) = (3.8) 

j=0 



Consequently for n = 0, 1, . . . , T — 1 

1 k 



i^(a^^+") = K{1) 



h{a 



T-l\ 



/i(a"-^). (3.9) 



Proof of theorem 3.1: Let t = kT + v, then it follows from (3.5) and (3.9) that 

E^ikT + V) - _ i^(l)[Ma^-^)]^Ma-+-^) ^ 
+ ^) - Kian) K(l)Ma"-^) ^ 

= [^(a^-^)]'=^(a^+"-^)[^(a"-^)]-^i?^(0) 
for A; = 0, 1, . . ., a > 1 and ra, f = 0, 1, . . . , T - !.□ 



Using theorem 3.1, one can easily verify that 

oo oo oo 

r=0 fe=0 fe=0 



= \h{a^+--')[h{a--')]-'\R^ (0) ^ |Ma^-^)|'= 

fc=o 

Now the process is short memory if h{a'^~^) < 1, otherwise the process is long memory. 



Remark 3.2 It follows from theorem 3.1 and relations (3.3) and (3.7) that R^{t), 



n = 0,l,...,T-l is fully specified by the values of {Rf{l), RfiO),j = 0, 1, . . . , T - 1}. 



By the following theorem, the necessary conditions are given to show that, the 2T elements 
{Rf{l), RfiO),j = 0, 1, ... ,r - 1} of the function -R^(r), given by (3.2) characterize the 
covariance function of the DT-SIM process. 

Theorem 3.3 The function R^ir) in relation (3.2) characterize the covariance function of 
a DT-SIM process, if these covariance functions satisfy the following condition 

{Rf ^ Rf mf+ii^) (3.10) 
/or j = 0, 1, . . . , T - 1 and R^ (0) = a'^'^"R^{0). 
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Proof: As by theorem 3.1 we have that the covariance function of the DT-SIM process is 
of the form (3.2), so for the proof of this theorem it is enough to show that every covari- 
ance function of the form (3.2) is the covariance function of a DT-SIM process. For this, we 
show in (i) that the covariance function R^{t) which is defined in (3.1) has scale invariance 
property with scale / = , i.e, satisfies the equality R^j^r^{kT + v) = f^R^ikT + v). We 
also prove in (ii) that R^iT) is the covariance function of a wide sense Markov process, as it 
satisfies in the relation (2.6). 

(i) According to (3.2) we have that 

R^^r{kT + v) = [^(a^-^)]^^(a^+-+^-i)[^(a"+^-^)]-^i?f+^(0) 
where by (3.3) and (3.7) 

„^,+n+T-l fif (1) T-rT-l fi.f (1) ™+n+r-l RfW 
/i(q!^+"+^-1) iij=0 Rf{0) i-i-j=0 RH(0)ilj=T ijH(o) 

h(an+T-l) " rin+T^l RfW " riT-1 RfW nn+T-1 Rfi~ 



and 



nv+n-l Rf+rW y-rv+n-l Rf(^) 
3=0 Rf^^m _ llj=o (0) _ /i(a''+'*-i) 



Thus R^+rikT + v)= P"R^{kT + v). 

(ii) By (3.1) and (3.2) we have that 

R^{kT + v)=Rx{a^^^''+'',a^) = [/t(a^-i)]*^/i(a''+"-^)[/i(a"-^)]-^i?f (0) 
= if(l)[^(«^-i)]^^(a''+"-i)[K(l)^(«"-i)]-ii2f(0). 

According to (3.8) 

R^{kT + v) = K{l)h{Q^'^+''+''-^)[K{l)h{a^~^)]-^R^,{d) 

= i^(a'=^+"+'')[i^(a")]-^i?^(0) = i^(a*^^+'^+")G(a") 
where by (3.4), (3.9) 

K{a^) = K{l)h{a^-% G{a^) = [K{a^)]-^R^ (0). 

Thus for T G W, a > 1 

RH^r) = Rx{a''+\a'') = K{a''+^)G{a''). 
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As a""*""^ > a"^, by Borisov condition (2.6) 

i^xK+^a") = G(mmK+^,a"))i^(max(a"+^a")) = 
Finally we show that G/K is positive and nondecreasing. Note that 

G(a") i?^(0) 



> 0. (3.11) 



Now we prove that F^a"^) = G{q^) / K{a^) under condition (3.10) is nondecreasing function. 
So we show that F(a^+i)/F(a'") ^1. Let r = A;T + n ^ and < n ^ T - 1, then we 
consider two cases, ^ n ^ T — 2 and n = T — 1. 

For ^ n ^ T - 2, by (3.11) and (3.9) we have that 



F(a'+i) _ iv-2(o«--") i?f^-+„+i(0) 



2 uH 



-Rn+l(0) 



F(«'-) i^2(afcT+n+l) iif^^jo) 

As X{-) is DT-SI with parameter space {a'^, A; G W} and scale ^ so by (3.1) 

RkT+n{^) = i?x(a'=^+",a'=^+") = £;[X(a*^^+'^)X(a'=^+")] 
= (a^)2'=^£;[X(a")X(a")] = a''^'^" {Q). 

Therefore 

i^(«^+^) ^ . /l(a"-^) 2 -^n+l(0) ^ ri^^(0)i2 -^n+l(0) ^ (0)^^^+! (0) 

F(a'-) ^ /i(a") ^ i?^(0) ^i?^(l)J i?^(0) (i?^(l))2 • 

Under condition (3.10) for j = n, F{a^~^^)/F{a'^) ^ 1 and F{a^) is nondecreasing. 
For n = T — 1 we have 

-Rf_i(o).2_R£(o)_ ^ i^-i(o)i^(o) 
^i?f_i(i)^ i?f_i(o) (i?f_i(i))^ ■ 

By a similar method and under condition (3.10) for j = T — 1, F{a^) is nondecreasing.D 
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4 T-dimensional self-similar Markov processes 



In this section we assume that {X{a"'),n G W} is a DT-SIM process with scale I = , and 
characterize the covariance function of the associated T-dimensional discrete time self-similar 
Markov process in theorem 4.1. 

Definition 4.1 The process U{t) = {U°{tso),U^{tsi),...,U^-'^{tsT-i)) for fixed positive 

SQ,Sl, 

• • • , st-1 and with parameter space T = {P; n G W, I > 1} is a T-dimensional discrete time 
self similar if the fallowings are satisfied 

(a) {U^ (t)} for all j = 0,1, ■■ ■ ,T — 1 is a DT-SS process with parameter space 

= {rsj;n G W}. 

(b) U^{-) and f/*(-) has self-similar correlation, that is 

Cov(C/^(Z"+"isj)C/*(Z"+"2si)) = l'^''" Cow {W{r'sj)U'{l''^ Si)) 

for every n, ni, n2 G W, i, j = 0, 1, • • • ,T — 1. 

Corresponding to the DT-SIM process {X{a^),n G W} with scale I = , a> 1 there exist 
a T-dimcnsional discrete time self-similar Markov process V{t) = (V^{t),V^{t), . . . , V^~^{t)) 
with parameter space T = {P; n G W, I > 1}, where I = , a > 1 

= := X(a"^+*^). (4.1) 

This is by the fact that F'^(P) are DT-SS process for A; = 0, 1, . . . , T - 1 according to the 
definition 2.5 for s = a'' and all n,ni,n2 G W 

Cov(W(r+'^i),F^(r+"2)) = Cov(X(a("+"i)^+^'),X(a("+"2)T+i)^ 

= a^'''^"CoY{X{a'''^+^),X{a'''^+')) = f''^Cov{V\l'''),V\V'^)). 

This equality provide the assertion (b) of definition 4.1 with Si = Sj = 1, i,j = Q, - ■ ■ , T — 1. 
Let Q^{n,T) = [Q^(n,T)]j_fe=o,i,...,r-i be the covariance matrix of then for j,k = 

0,1,...,T-1 

Qg(n,r) = E[V\V+^)V^{V)] = £;[X(a("+^)^+^)X(a'*^+*^)], r G W. 

Theorem 4.1 Let {X{a^),n G W} be a DT-SIM process with the covariance function {t) 
and let {T^(Z"),n G W}, defined in (4.1), be its associated T-dimensional discrete time self- 
similar process with covariance function Q^{n, r). Then 

Q'^in, t) = a^'^''^Q"{T) = a2«^^C^^i?^[/i(a^-^)]^ r G W (4.2) 
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where h{-) is defined by (3.3) and the matrices Ch and Rh are given by 

^O.T-l 



Ch 



'^00 '^01 



^10 
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Cfk = ^(a^-^)[/i(a*=-^)]-S and 



Rh = 



<(0) 
i?f(0) 








••• i?T_i(0) 

Proof: As X{-) is DT-SI with parameter space {a^ , k G W} and scale , so by (3.1) 

Q%{n,r) = a'-^^E[X{a^^+^)X{a')] = a'-^^R^{rT + j - k). 
Let i?f (rT + j - k) = Q%{t), then 

Q%{n,r) = a--^^Q%{r) 

hence 



R^{tT) 
R^irT + l) 



R^{tT - 1) 
R^irT) 



R^_,{tT-{T-1)) 
R^_,{tT-{T-2)) 



RI}{tT + T -I) R^{tT + T-2) 

and by the Markov property of X(-) from theorem 3.1, for O^j — A;^r — 1 
R^irT + j -k) = \Ji{a^-')Yh{a^-')[h{a^-')]-'R^{G). 

jk 



Let = h{a^-^)[h{a''-^)]-\ so 



QS(r) = [Ma^-^)]-Cf,<(0) 



T-l\-\Tr^H r>H, 



T G W. 



Thus we can represent Q^{n,T) as 



T-lMT 



Coo^o^(O) 



Co^i?f(0) 
CdR^O) 



Cnrp_^ -Rt^_1 (0) 



^O^T-l-"'r- 

'^i,r-i-^r-i(0) 



Ct-i,o^o{^) C'^-i.i-^f (0) ■■■ C'T-i,r-i-^T-i(0) 
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